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Abstract. The propagation of a modulated beam of charged particles in a planar 
crystal channel is investigated. It is demonstrated that the beam preserves its 
modulation at sufRciently large penetration depths which opens the prospect of using 
a crystalline undulator as a coherent source of hard X rays. This finding is a crucial 
milestone in developing a new type of lasers radiating in the hard X ray and gamma 
ray range. 
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In this Letter we study for the first time the evolution of a modulated particle beam 
in a bent planar crystal channel and demonstrate that it preserves its modulation at 
sufficiently large penetration depths which opens the prospect of using a crystalline 
undulator as a coherent source of hard X rays. Solving this problem is of crucial 
importance in the theory of the Crystal Undulator based Laser (CUL) [H [2], [3] — a 
new electromagnetic radiation source in hard X and gamma ray range. 

Channeling takes place if charged particles enter a single crystal at small angle 
with respect to crystallographic planes or axes [4j. The particles get confined by the 
interplanar or axial potential and follow the shape of the corresponding planes and axes. 

A single crystal with periodically bent crystallographic planes can force channeling 
particles to move along nearly sinusoidal trajectories and radiate in hard X and 
gamma ray frequency range (see Fig. [1]). The feasibility of such a device, known as 
the 'crystalline undulator', was demonstrated theoretically a decade ago [1] (further 
developments as well as historical references are reviewed in [5]). More recently, an 
electron based crystalline undulator has been proposed [6]. 




Figure 1. Schematic representation of the crystalUne undulator. 

The advantage of the crystalline undulator is in extremely strong electrostatic fields 
inside a crystal which are able to steer the particles much more effectively than even 
the most advanced superconductive magnets. This fact allows to make the period 
of the crystalline undulator in the range of hundreds or tens micron which is two to 
three orders of magnitude smaller than that of conventional undulator. Therefore the 
wavelength of the produced radiation A ~ Au/(27^) (7 ~ 10^-10^ being the Lorentz 
factor of the particle) can reach the (sub)picometer range, where conventional sources 
with comparable intensity are unavailable [7]. 

Even more powerful and coherent radiation will be emitted if the probability density 
of the particles in the beam is modulated in the longitudinal direction with the period 
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A, equal to the wavelength of the emitted radiation. In this case, the electromagnetic 
waves emitted in the forward direction by different particles have approximately the 
same phase [8]. Therefore, the intensity of the radiation becomes proportional to the 
beam density squared (in contrast to the linear proportionality for an unmodulated 
beam). This increases the photon flux by orders of magnitude relative to the radiation 
of unmodulated beam of the same density. The radiation of a modulated beam in an 
undulator is a keystone of the physics of free-electron lasers [9j. It can be considered 
as a classical counterpart of the stimulated emission in quantum physics. Therefore, if 
similar phenomenon takes place in a crystalline undulator, it can be referred to as the 
lasmg regime of the crystalline undulator. 

The feasibility of CUL radiating in hard X ray and gamma ray range was considered 
for the fist time in [H |2] . Recently, a two-crystal scheme, the gamma klystron, has been 
proposed [3]. 

A simplified model used in the cited papers assumed that all particle trajectories 
follow exactly the shape of the bent channel. In reality, however, the particle moving 
along the channel also oscillates in the transverse direction with respect to the channel 
axis (see the shape of the trajectory in Fig. [1]). Different particles have different 
amplitudes of the oscillations inside the channel (Fig. [2], upper panel). Similarly, the 






Figure 2. Due to different amplitudes of channeling oscillation (upper panel) and 
different momentum directions in the (xz) plane (lower panel) , the initially modulated 
beam gets demodulated. Open and filled circles denotes the same particles at the 
crystal entrance and after travelling some distance in the crystal channel, respectively. 



Stable propagation of a modulated particle beam in a bent crystal channel 



4 



directions of particle momenta in (xz) plane are slightly different (Fig. [21 lower panel). 
Even if the speed of the particles along their trajectories is the same, the particles 
oscillating with different amplitudes or the particles with different trajectory slopes 
with respect to z axis have slightly different components of their velocities along the 
channel. As a result, the beam gets demodulated. An additional contribution to the 
beam demodulation comes from incoherent collisions of the channeling particles with 
the crystal constituents. 

In the case of an unmodulated beam, the length of the crystalline undulator 
and, consequently, the maximum accessible intensity of the radiation is limited by the 
dechanneling process. The channeling particle gradually gains the energy of transverse 
oscillation due to collisions with crystal constituents. At some point this energy exceeds 
the maximum value of the interplanar potential and the particle leaves the channel. The 
average penetration length at which this happens is known as the dechanneling length. 
The dechanneled particle does not follow the sinusoidal shape of the channel any more 
and, therefore, does not contribute to the undulator radiation. Hence, the reasonable 
length of the crystalline undulator is limited to a few dechanneling lengths. A longer 
crystal would attenuate rather then produce the radiation. Since the intensity of the 
undulator radiation is proportional to the undulator length squared, the dechanneling 
length and the attenuation length are the main restricting factors that have to be taken 
into account when the radiation output is calculated. 

In contrast, not only the shape of the trajectory but also the particles positions with 
respect to each other along z axis are important for the lasing regime. If these positions 
become random because of the beam demodulation, the intensity of the radiation 
drops even if the particles are still in the channeling mode. Hence, it is the beam 
demodulation rather than dechanneling that restricts the intensity of the radiation of 
CUL. Understanding this process and estimating the characteristic length at which this 
phenomenon takes place is, therefore, a cornerstone of the theory of this new radiation 
source. 

Let us consider the distribution f(t, s; ^, Ey) of the beam particles with respect to 
the angle between the particle trajectory and axis z in the {xz) plane ^ = arcsinp^/p ^ 
Px/p and the energy of the channeling oscillation Ey = p^y/2E + U{y) [||- Here p, px 
and Py are, respectively, the particle momentum and its x and y components, U (y) is 
the interplanar potential which we approximate by a parabola U{y) = ?7max(y/ymax)^, 
and E is the particle energy (we will consider only ultrarelativistic particles, therefore 
E ^ p). It can be shown that the evolution of this distribution in the crystal channel 
with the time t and the longitudinal coordinate (penetration depth into the crystal 
along the curved channel) s can be described by the following differential equation of 
Fokker- Planck type: 



I We chose the system of units in such a way that the speed of hght is equal to unity. Therefore, mass, 
energy and momentum have the same dimensionahty. This is also true for length and time. 




(1) 
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Here Dq is the diffusion coefficient that is dominated by the scattering of the beam 
particles by lattice electrons. The particle longitudinal velocity averaged over the 
undulator period, Vs, is given by 

'l---^-^^ (2) 
272 2 2^7 ■ ^ ^ 

Equation ([T]) is akin to the equation describing dechanneling process (see e.g. [lU]) and 
is derived in a similar way. The novel feature of it is the presence of time variable, 
which allows to describe time dependent (modulated) beams. Additionaly, it takes into 
account scattering in the {x, z) plane. 

If the beam is periodically modulated (bunched) the distribution can be represented 

Fourier series: 



00 



f{t,s;^,Ey) = gj{s;^,Ey)exp{ijujt). (3) 



j=-oo 



with g*{s]^,Ey) = g^j{s] C,, Ey) to ensure the real value of the particle distribution. 
Since Eq. ([1]) is linear, it is sufficient to consider only one harmonic. Substituting 

f{t, s; ^, Ey) = g{s; ^, Ey) exp{iujt) one obtains 



dg 

iujg{s; ^, Ey) + —Vs = Dq 



d f dg \ 1 d g 



dEy y^'^oEyj ' Ede\' ^'^^ 

To simplify this equation, we make the substitution g{s; ^, Ey) = exp {—ius) g{s; ^, Ey) 
and assume that the variation of g{s; ^, Ey) within the modulation period is small: 
dg/ds -C ujg{s; ^, Ey). This allows us to neglect the terms {l—Vs)dg/ds while keeping the 
terms (1 — Vs)ojg{s] ^, Ey). The resultant partial differential equation for g[s; ^, Ey) can 
be solved by the method of separation of variables. Putting g{s] ^, Ey) = S{s)E{^)S{Ey), 
we obtain a set of ordinary differential equations: 

Do 1 d'Ejo . e . 

=Cc, (5) 

D, d [ dS{Ey) \ Ey _ 



S(^Ey) dEy y dEy J 2E 
' '^^'K^ =Cg, (7) 



S{s) ds 272 

where Cs, Q and Cy do not depend on any of the variables s, C, and Ey and satisfy the 
condition 

= Q + Cy. (8) 

Eq. has the form of the Schrodinger equation for the harmonic oscillator. Its 
eigenvalues and eigenfunctions are, respectively. 



Q,. = -(l + ^)y^(^n+^), n = 0,l,2,... (9) 

and 



2D, 








IniO = Hn I e 1 exp I -^JZ^e 1 • (10) 
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Here -ff„(. . .) are Hermite polynomials. 

Eq. ([6]) can be reduced to the Laguerre differential equation, so that its solution 
can be represented as 



S.(E,)^eM-'-i^Jj^E,\ Ud+id^E,) (11) 



where Ly{. . .) is the Laguerre functioii§| and Uk is related to the eigenvalue Cy^k via 

Cy,k = -^-^\l^{2vk + l), fc = 1,2,3,... (12) 

The eigenvalues can be found by imposing the boundary conditions. The maximum 
energy of channeling oscillations in a bent channel with parabolic potential is given by 
^max _ Uraax{^ — CY , whcre C = -f^c/^max i^ defined as a ratio of the centrifugal force Fc 
to the maximum value of the interplanar force t^^ax 0- Therefore, the density boundary 
condition has the form: 



L., [{l+^)^^U^^^{l-Crj=0. (13) 

Equation (fT3|) has to be solved for u^. (the subscript k enumerates different roots of the 
equation) and the result has to be substituted into (|T2l) . 

It is convenient to represent the eigenvalues in the form 

Cy,, = -^J^^-tu;elP,iK,C). (14) 

Here = 4f/max/(jo,i-^o) dechanneling length in a straight channel [10] (jo,fc is 

fc-th zero of the Bessel function Jo{e)), = ^y2U^ax/E is the corresponding Lindhard's 
angle. We introduce the parameter 

= ir^el (15) 

where A = 2tt/uj is the spatial period of the modulation. The functions ak{K,,C) and 
Pk{i^,C) can be represented as 

ak (/t(l-C )^) 

[i-cy 

Pk{f^,c) = ii-cYPk{<i-cY) (17) 

Here akin) = a^^n^C = 0) and Pk{f^) = Pk{i^,C = 0) are the cooresponding functions 
for the straight channel which are found by solving numerically Eq. (JT3l) combined with 

m- 

Using (IH]) , one finds the solution of Eq. ([7]) : 



ak{fi,C) = — — (16) 



d 



— ZUJS 



Jo,i, 

27 2jo,iV/« 



§ At nonnegative integer values of the Laguerre function is reduced to the weU known Laguerre 
polynomials. In the general case that is relevant to our consideration, it can be represented by an 
infinite series: L^(e) = Y^TLo Umloi"^ ~ i')eV(jO^- 
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Hence, the solution of Eq. (jl]) is represented as 

oo oo 

g{s] ^, Ey) = exp {-iujs) Cn,k^n{0£kiEy)Sn,k{s), (19) 



n=0 k=l 



where the coefficients Cn,k are found from the particle distribution at the entrance of 
the crystal channel. Due to the exponential decrease of Sn,k{s) with s (see ( ITSl) ). the 
asymptotic behaviour of g{s]^,Ey) at large s is dominated by the term with n = 
and k = 1 having the smallest value of the factor [ak{K,,C) + {2n + 1)y^/jo,i] in the 
exponential. Therefore, at sufficiently large penetration depths, the particle distribution 
depends on s as g{s] ^, Ey) oc exp (— s/Ldm — iuj/us s) where L^m is the newly introduced 
parameter — the demodulation length: 



and Us is the phase velocity of the modulated beam along the crystal channel 

1 



(20) 



(21) 



This parameter is important for establishing the resonance conditions between the 
undulator parameters and the radiation wavelength. It will be analysed elsewhere. 




Figure 3. The ratio of tlie demodulation length Ldm (f20| to the dechanneling length 
in the straight channel versus the parameter k (fTS]) for different values of Tsyganov's 
factor C. The corresponding asymptotic values at k — > are shown by thin horizontal 
lines. 

In this Letter we concentrate our attention on the demodulation length. This 
parameter represents the characteristic scale of the penetration depth at which a beam 
of channeling particles gets demodulated. 

Fig. [3] presents the dependence of the ratio L^^/ on the parameter k. At k — 0, 
the demodulation length approaches (1 — C^L^ which is the dechanneling length in the 
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bent crystal. It is seen that the demodulation length is smaller than dechanneling length 
by only 20-30% at k < 1 for C varying between and 0.3. This is an important result, 
since it was proven for a number crystals channels [5J that the dechanneling length 
is sufficiently large to make the crystalline undulator feasible. According our present 
result, the demodulation process does not put essential additional restrictions on the 
undulator length and, therefore, the crystalline undulator can generate also coherent 
radiation provided that it is fed by a modulated beam. On the contrary, the ratio 
noticeably drops for k > 10. Therefore, further investigations are needed to clarify 
whether the lasing regime is feasible in this k range. 

The dependence of the parameter n on the energy of the emitted photons, hLj = 
27rh/X, is shown in Fig. |H The calculation was done for 1 Gev[Jj] positrons using the 
formula for the dechanneling length from [Sj [10] . 
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Figure 4. The parameter k (fT5)) versus the photon energy huj for different crystals 
and crystallographic planes. 

As one sees from the figure, there exist such crystal channels that k, = 1 corresponds 
to hw = 100 — 300 keV. Therefore, CUL is most suitable for the application in this 
photon energy range. The lower limit on the photon energy is dominated by a very 
strong absorption of the radiation in the crystal, while the upper limit on the photon 
energy is set by the decreasing value of Ldm- 

According to our estimations, the brilliance as high as 
]^g25_]^g26 Photous/ (s mm^ mrad2 0.1%BW) can be obtained in a CUL fed by a com- 
pletely modulated positron beam with current 1 kA and particle density 10^^ cm~^. 

In conclusion, we have studied the propagation of a modulated particle beam in 
a bent planar crystal channel. It has been demonstrated that the beam preserves 

II Note that k depends weakly (logarithmically) on the particle energy. Therefore, changing the beam 
energy by an order of magnitude would leave Fig. |4] practically unchanged. 
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its modulation at the penetration depths which are sufficient for producing coherent 
radiation with the photon energy of hundreds keV. This opens the prospects for creating 
intense monochromatic radiation sources in the frequency range which is unattainable 
for conventional free electron lasers. Developing suitable methods of beam modulation 
would be the next milestone on the way towards this goal. 
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